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ABSTRACT 

We present a global analysis of kinematics and metallicity in the nearest SO galaxy, NGC 3115, 
along with implications for its assembly history. The data include high-quality wide-field imaging 
from Suprime-Cam on the Subaru telescope, and multi-slit spectra of the field stars and globular 
clusters (GCs) obtained using Keck-DEIMOS/LRIS and Magellan-IMACS. Within two effective radii, 
the bulge (as traced by the stars and metal-rich GCs) is flattened and rotates rapidly (v/a > 1.5). At 
larger radii, the rotation declines dramatically to v/a ~ 0.7, but remains well-aligned with the inner 
regions. The radial decrease in characteristic metallicity of both the metal-rich and metal-poor GC 
subpopulations produces strong gradients with power law slopes of —0.17 ±0.04 and —0.38 ±0.06 dex 
per dex, respectively. We argue that this pattern is not naturally explained by a binary major merger, 
but instead by a two-phase assembly process where the inner regions have formed in an early violent, 
dissipative phase, followed by the protracted growth of the outer parts via minor mergers with typical 
mass ratios of ~ 15-20:1. 

Subject headings: galaxies: bulges — galaxies: formation — galaxies: halos — galaxies: individual: 
NGC 3115 — galaxies: kinematics and dynamics — globular clusters: general 



1. INTRODUCTION 

The faint outer regions of galaxies are particularly use- 
ful probes of the earliest phases of galaxy formation. 
Here relaxation times are long, and subtle formational 
clues are preserved not only in their integrated stellar 
light but also in bright tracer populations such as glob- 
ular clusters (GCs) and planetary nebulae (PNe). 

Recently, wi de field observations of integrated stellar 
light and PNe (jProctor et al.l 120091 : iCoccato et aTll2009[ ) 
have revealed that the relative kinematic homogeneity 
among the central reg ions of early-type galaxies (ETGs) 
(Emsellem et al. 2007) gives way to surprising diversity 
in the outlying parts, e.g., central fast rotators can some- 
times rotate slowly at larger radii. 

The implication is that inner and outer regions 
are somewhat decoupled. In the classic gas-rich bi- 
nary major- merger scenario for forming galaxy bulges, 
the spatial dependence of dissipational processes could 
produce dramatic k inematical transitions in radius 
([Hoffman et~a!1l20lol ). 

Alternatively, inner and outer spheroids may form in 
two separate phases, which could explain the strong size 
and intern al density evolution of ETGs from z y 2 t o 
today (e.g. iDaddi et al.l 120051 : Ivan Dokkum et al.l 120081 ). 
This inside-out growth hypothesis means that the inner 
regions form early as compact stellar spheroids while the 
outer parts grow later through the accretion of smaller 
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galaxies (e.g. JNaab et al.l 120091 : fOser et al.ll201QD . 

One way to distinguish between the merger and inside- 
out scenarios is to search for their signatures in the 
properties of GC sy stems. In their archetypal study, 
iSearle fc Zinnl (jl978[ ) used the joint distribution of spa- 
tial, kinematical, age, and metallicity properties of GCs 
to infer the hierarchical buildup of the Milky Way's outer 
halo. A key concept here is that building galaxies via the 
stochastic infall of satellites leads to flattened metallicity 
gradients at large radius, an e ffect now observed within 
the G C syst ems of two ETGs ([Harrisl l2009l: iForbes et all 
120111 : cf also lCoccato et al.ll2010f ). Likewise, in wet major 
merger remnants gas dissipation should create metallic- 
ity gradients in the central regions, with flatter profiles 
in the outer regions where an y initial grad i ents have been 
mixe d up in the merger (e.g., [White 1980; Hop kins et all 
2009). The similar expectations for large-scale metal- 
licity structure from each of these scenarios necessitates 
additional information, such as from kinematics, to place 
such results in the proper context. 

A useful complication with the GCs is their bimodality, 
as each galaxy generally has a metal-poor and a metal- 
rich subpopulation (MPGCs and MRGCs). In the Milky 
Way these are associated with the stellar halo and the 
bulge or thick disk, respectively, while in ETGs there are 
stronger MRGC components reflecting their more dom- 
inant bulges. These subpopulations permit the investi- 
gation of two distinct phases of galaxy assembly (stellar 
halo and bulge) far beyond the Local Group. 

As part of our ongoing surveys of ETGs we have stud- 
ied the nearest giant SO galaxy, NGC 3115 (distance 9 
Mpc, effective radius R e =bl"=2.b k pc, B/T > 0.9, in- 
clination 86°: iCapaccioli et al.l Il987h . in unprecedented 
photometric and spectroscopic detail. In particular, we 
have obtained extensive kinematic data-sets of its field 
stars and GCs. Here we map out the global kine- 
matic and metallicity structure of its dominant, extended 
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Fig. 1. — Spatial and kinematic maps of NGC 3115. Top: i-band image (a) with SKiMS (stellar light) measurement locations (purple 
stars). Maps of stellar rotation (b), velocity dispersion (c) and /13 (d), each sharing the same spatial scale as (a). Bottom: Zoomed-out 
view (e) with spectroscopically confirmed metal-rich (red circles) and metal-poor GC (blue crosses) locations. The inset rectangle denotes 
the region shown in (a) where selected i-band isophotes are also overplotted. Rotation maps, made using the MRGC (f) and MPGC 
(g) datasets, with the same spatial scale as (e), and overlaid by an outline of the zoomed-in stellar fields from (b)-(d). The kinematic 
maps were created using kriging techniques for statistical ly reconstructing sparsely-sampled functions, using "interpolation" and "fitting" 
approaches for the stellar and GC datasets, respectively (Press et al. 2007). The minimum/maximum color-bar scaling is recorded in the 
bottom right corner of each panel (all in units of km s — 1 except for (d), which is dimensionless). While these maps effectively convey 
the broad-brush kinematic behavior, we note that interpolating sparse data sets may lead to over-smoothing and anomalous small-scale 
structure, particularly at larger radius. 



bulge, and examine some implications for its assembly 
history. 

In $2] we describe our imaging and spectroscopic obser- 
vations. Kinematic and metallicity profiles are presented 
in 331 We discuss how these results constrain the forma- 
tion mechanisms of the outer bulge and halo in § [4j and 
summarize the conclusions in § [5l 

2. OBSERVATIONS 

Here we briefly summarize the observations used in 
the present analysis; complete details will be presented 
in a follow-up paper. We obtained gri images with 0.5"- 
0.7" seeing using Suprime-Cam on the 8.2 m Subaru tele- 
scope. Sections of this imaging are shown in Fig. la,e to 
provide orientation for the spectroscopic measurements 
discussed below. The photometry allows us to select GC 
candidates for follow-up spectroscopy and to character- 
ize the distributions of number density and metallicity 
for the GC system. Color is used as a proxy for metal- 
licity since GCs are generally old (see Brodie fc Straderl 
2006; iKuntschner et all 12002). The color boundary we 
use for MPGCs/MRGCs in NGC 3115 is (g - i) = 0.91 
mag ([Fe/H] = -0.8 dex). 

We performed spectroscopy using DEIMOS on the 
Keck-II telescope to focus o n the calcium triplet re- 
gion (~8580 A, R ~ 5200; see iRomanowskv et all [20091 : 
iFoster et al.ll20lil ). We have confirmed 150 GCs and also 
obtained spectra of the background galaxy light in 166 
slits for use in a novel technique ca lled "SKiMS" (Stel- 
lar Kinematics fr om Multiple Slits: iProctor et al.l 2009; 
IFoster et al.ll20lH ). Another 15 and 11 unique GC spec- 
tra were obtained using LRIS (-3500-7500 A, R ~ 1100- 
2200) and IMACS (-4000-7000 A, R ~ 2200), on the 
Keck-I and Magellan telescopes, respectively. 

GC velocities were measured by cross-correlating their 
absorption-line spectra with stellar templates. We im- 
posed a lower velocity cut at 350 km s _1 to remove 
foreground stars. Fainter SKiMS spectra were coad- 



ded to reach a tar get S/N of 25 using the V oronoi 2D- 
binning method by Cappe llari fc Copinl (|2003l ). The stel- 
lar line-of-sight velocity distribution of each spectrum 
was parameterized as a tr uncated Gauss- Hermite series 
(V, cr, ft 3 , ft 4 ) using pPXF dCappellari fc Emsellemll200l 
and 13 template stars, with parameter uncertainties es- 
timated using Monte Carlo modeling of data-sets with 
added noise. 

We supplement our SKiM S catalog w it h lon g-slit 
stellar kinematics data from iNorris et al.l ([2006) and 
llllingworth fc Schechterl (|1982l ). and incorporate 20 GC 
velocities from pr evious studies (jKuntschner et al.l 120021 : 
iPuzia et al J 12004( 1 

3. RESULTS 

Fig. 1 shows two-dimensional maps of the kinematic 
data. In the top row, the thin disk is evident as a narrow 
stripe within an extended and rapidly rotating bulge. 
The thin disk's lowered velocity dispersion, and anti- 
correlated V and /13, extend well into the bulge, suggest- 
ing an additional embedded component like a very thick 
disk. On larger scales (Fig. 1 bottom row), distinct rota- 
tion is seen in the MRGC and MPGC subsystems, while 
there is a general trend for the rotation to weaken in the 
outer regions. 

To simplify the rotation and dispersion trends into 
one-dimensional profiles, we use a variation of the 
"kin emetry" technique dev eloped for the SAURON sur- 
vey (jKrajnovic et al.l 120061 ) as optimized for data with 
discrete spatial and v elocity sampling (jProctor et al.l 
120091 IFoster et alJl20lJ l. This method samples the kine- 
matic field (e.g., rotation) using concentric elliptical an- 
nuli, and fits the data to flattened sinusoidal models as 
a function of position angle. For the discrete velocity 
data (GCs), rotation and dispersion are fitted simulta- 
neously through a maximum likelihood method (Figs. 
2a-d). The position angles and ellipticities of the rota- 
tion field and the sampling bins (PAkin, tkin) are part of 
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Fig. 2. — Rotation models for NGC 3115. Panels (a)-(d) plot line-of-sight velocities in two radial bins (R/R e ~ 1, 5) versus position 
angle for MRGCs (top) and MPGCs (bottom). Solutions for rotational velocity (V TO t: solid black curves) and velocity dispersion (a) are 
determined for each dataset using the kinemetry method described in the text (dotted curves denote V vo t ± 2a). Smoothed rolling- fit 
rotational profiles with uncertainty envelopes are shown for the stellar light (black curve) and MRGCs (red curve) in (e), and for the 
MPGCs (blue curve) in (f). In (g) the v/ a profile of NGC 3115' s MRGC subpopulation (red curve) is compared against a simulated merger 
remnant with a 1:10 mass-ratio (dashed: Bournaud et al. 2005). (h) contains the z-band ellipticity profile. 



the fit for the SKiMS data but are not well constrained 
for the discrete velocity data, which we assume follows 
the stellar isophotes (PA=43.5°, e=0.5). Our results are 
insensitive to reasonable variations in these parameters. 
Uncertainties are estimated via Monte Carlo fitting of 
mock data-sets. 

The resulting rotation profiles for the different sub- 
components are shown in Figs. 2e and 2f, where rolling 
fits with radius are used to capture the details of any 
radial kinematic transitions (e.g from the inner to outer 
bulge/halo). Within ~ 1.5 R e the MRGC system rotates 
nearly as rapidly as the stellar bulge, supporting the co- 
evolution of these two components, as also inferred from 
their similar ages and metallicities (jNorris et al.l 12006). 
At larger radii, this rotation decreases dramatically (see 
also Figs. 2a, b). The MPGCs have moderate rotation 
with a decline outside ~ 4i? e . 

An alternative rotation profile for the MRGCs is shown 
in Fig. 2g, after normalizing by the local velocity disper- 
sion. The photometric ellipticity profile is also plotted, 
showing a decrease with radius that parallels the rota- 
tional gradient. The overall implication is for a bulge 
that has a high degree of rotational flattening in its 
central regions, while becoming rounder and dispersion- 
dominated in its outskirts. 

Having found kinematic transitions in both GC sub- 
populations, we look for analogous transitions in the ra- 
dial metallicity profiles. First we summarize the overall 
color distribution of the GCs in Fig. 3a, which shows a 
classic bimodality. We will assume that this bimodality 
persists with increasing radius, but that the location of 
the color peaks may shift. At large radii we must cope 
with the contaminating effects of foreground stars, whose 
color distribution we also show in Fig. 3a, and use to con- 
struct Monte Carlo mock datasets to iteratively correct 
for the contaminant bias on the color peak locations. 

Fig. 3c shows color versus radius, both for individ- 
ual GC candidates and for the fitted peak locations. 
Both GC subpopulations have radially-decreasing col- 
ors, which we quantify as power-law color gradients with 



slopes of —0.05 and —0.07 mag per dex for MPGCs and 
MRGCs, respectively. Using our own empiri cal calibra- 
tion t o the (g — z) color used in ACS surveys (jPeng et al.l 
2006), the gradients are —0.07 an d —0.10 mag per dex. 
Converting to [Fe/H] metallicity (Pe ng et al.l l2QQ6h . we 
estimate gradients of -0.38 ± 0.06 and -0.17 ± 0.04 dex 
per dex. 

To our knowledge, this is the first time that metal- 
licity gradients in both GC subpopulations have been 
measured to large radii in any galaxy besides a few very 
massive ellipticals (see Section [1]). It is also one of the 
first cases of any galaxy type where joint rotation and 
metallicity gradients are observed in t he halo (see also 
NGC 4697: iMendez et all I2005L [2009; and NGC 4125: 
IPu et al.ll20T(l . 

4. IMPLICATIONS 

We now consider some possible implications of the ro- 
tation and metallicity gradients for NGC 3115's assem- 
bly history. The central bulge properties are generally 
consistent with a standard major merger picture, with 
the very high amount of rotation in this case indica- 
tive of a gas-rich merger with an uneven mass-ratio (e.g., 
iBournaud etHI 120051 : iNaab et al.ll2QQ6D . Alternatively, 
the inner bulge might have formed via the inward migra- 
tion of giant star forming clumps within a turbulent disk 
fed b y cold streams f rom the cosmic web at early epochs 
(e^ lNoguchil 119991 : lElmegreen et al.l 120081 : iDekel et al.l 
2009). In either case, the exceptionally high inner-bulge 
rotation in NG C 3115 may require a residual thick dis k 
component (cf. iCretton et all 120011 : iWuyts et al1l201Q[ ). 
Also, the s urvival of an old th in disk with a mean age of 
~ 5-8 Gyr (jNorris et a l. 2006) means that most of the ac- 
tion in the central regions occurred before z ~ 0.5-1, and 
that there have been no major mergers more recently. 

At larger radii there are generic expectations from ma- 
jor mergers for rotation profile behavior. These rem- 
nants are generally expected to have rapid outer rota- 
tion resulting from both residual disk spin and the con- 
version of orbital into internal angular momentum (e.g., 
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Fig. 3. — Color distributions (g — i)o of point sources around NGC 3115. Left: Overall photometrically-selected GCs (open histogram) 
and of stellar contaminants (filled; as measured at R/R e > 14). Middle: Resulting true GCs (open) after contaminant subtraction. Also 
shown are the spectroscopically confirmed GCs (filled). Right: Color (left-axis) and [Fe/H] metallicity (right-axis; see text for details) 
versus radius. Individual GCs are shown as squares according to the legend at right. Solid curves denote r olling fits of the peak color of 
each GC subpopulation as estimated from a combination of mixture modeling methods ( Ashman et al. 1994) and an iterative Monte Carlo 
scheme to correct for the contaminant color-bias. Associated error envelopes reflect both statistical and photometric uncertainties. We 
convert the peak color difference between the inner (1-2 R e ) and outer (> 8 R e ) portions of each GC system into merger mass-ratios (see 
text), with sample values shown on the right-hand side. 



iBendo fc Barnesll2QQrt ICretton et alJl2Q01h . in stark con- 
trast to the declining rotation observed in NGC 3115. To 
make this difference more explicit, we searched through 
the literature for extended rotation profiles from simu- 
lations of major merger remnants, choosing a represen- 
tative example that comes close to reproducing the cen- 
tral rotation of NGC 3115's bulge and MRGC system. 
The chosen remnant is the result o f a 1:10 spiral-spiral 
merger from Bou rnaud et al.l (|2005l ). which we overplot 
in Fig. 2g, showing the discrepancy in the outer regions 
between the high rotation predicted and the low rotation 
observed. 

Gas-rich 1:1 merger remnants with small percen- 
ters can also produce declining outer rotation, but 
the outer "dry" part of the remnant is generally ex- 
pected to sh ow kinematical misa lignment with the in- 
ner regions (Hof fman et al.l I2Q10L and in preparation). 
Possible example s of this scenario include NGC 5128 
(jPeng et alJl200l and NGC 4125 (jPu et aLll2QToh . but 
other cases like NGC 3115 with decreasing but well- 
aligned rotation suggest there must be another explana- 
tion (cf. NGC 821 a nd NGC 3377: iProctor et aLl[200l 
iCoccato et al.l l2QQ9b and earl y arguments along these 
lines bv lScorza fc Benderlll995l ). 

Without exhaustive simulations of major mergers, we 
cannot rule out the possibility that finely tuned param- 
eters (viewing angle, impact parameter, etc.) would 
reproduce the observed kinematics of such systems. 
Nonetheless, it seems more natural to consider a two- 
phase assembly scenario (Section [T]) in which ETGs form 
inside-out. In this case, inner bulges form at high redshift 
while subsequent outer bulge and halo growth is driven 
primarily by dry minor- merger accretion events. The 
satellites fall in from many differe nt directions and pro- 
vide little net rot ational support (IVitvitska et al.l [ 2002: 
lAbadi et al.1l2QQ6b iBournaud etHll2QQ7t IQu et al.ll2Q10h . 



The radial decline in rotation of the MRGC system could 
then represent a transition from an inner bulge formed in 
violent, dissipative processes at high redshift, to an outer 
spheroid (around one third of the bulge mass in the case 
of NGC 3115) built largely from accreted material over 
a more protracted period. 

The MPGCs also show a marked rotational decrease, 
albeit at a larger radius. Here the theoretical picture is 
less clear, but we postulate an inner metal-poor stellar 
halo formed in-situ at high z followed by the accretion 
of outer material t hat also creates the outer bulge (cf. 
IZolotov et aTll2009[ ). 

Declining metallicity profiles are expected in this two- 
phase assembly scenario since the lower-mass accreted 
systems should be more metal-poor than the central 
galaxy (jNaab et al.l 120091 : iBezanson etall I2009D . The 
prediction (see JQ) is of a downward transition from an 
inner metallicity profile, whose shape reflects the detailed 
in-situ star formation history, to an outer profile flattened 
by radial mixing and primarily composed of accreted ma- 
terial. We assume this holds for both the stellar and GC 
metallicity profiles, though systematic offsets may arise 
from the difference between the ensemble properties of 
stars and GCs within incoming satellites (e.g. GCs are 
typically at larger galactocentric radii). As there are not 
yet any quantitative predictions, we adopt a schematic 
model where GC peak metallicities (or colors) are mark- 
ers of their host galaxies' masses (or luminosities), using 
the known correlations between these parameters at low 
z (see below). The peak GC color in the outer bulge/halo 
then indicates the characteristic luminosity of the ac- 
creted systems. 

The color-mass relations may well evolve with time, so 
we take a more general approach of considering the GC 
color difference between the central and outlying regions 
as an indicator of the characteristic mass-ratio that as- 
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semb led the outer gala xy. Using the color-mass relations 
from Peng et~aH (|2QQ6h for the MPGCs and MRGCs sep- 
arately, we then generate predicted outer GC peak colors 
for various stellar mass-ratios \\x in NGC 3115, overplot- 
ting these in Fig. 3c. For the MPGCs we find x = 200 
(15 < x < 2900; the color-mass relation is shallow so this 
constraint is weak) and for the MRGCs we find x = 8 
(3< x <19). In combination, we estimate x ~ 15-20. 

Very roughly, the accreted galaxies had luminosities of 
^ 5 x 1O 8 L^ )0 , equivalent to a dwarf elliptical. Such 
galaxies tod ay each host ~ 20 MPGCs and - 3 MRGCs 
on average (jPeng et al.l [2008). and could therefore ac- 
count for a total of - 350 MPGCs and - 50 MRGCs in 
the outer regions of NGC 3115. These numbers are differ- 
ent from the ~ 260 and ~140 GCs that we estimate NGC 
3115 to have outside 2 R e , but this discrepancy is not 
significant given the uncertainties in these calculations. 
A fundam ental complication arise s from biased galaxy 
assembly (jBrodie fc StradeH 12006). where the low- mass 
galaxies that were accreted at high z are thought to have 
hosted GCs of higher metallicity than their present-day 
counterparts in lower density environments. This means 
that the galaxies accreted by NGC 3115 may have had 
somewhat lower masses than in our simplified calcula- 
tion. 

As a demonstration of the ongoing accretion process, 
NGC 3115 does have a dwarf companion of this lumi- 
nosity at 45 kpc project ed distance, NGC 3115B. This 
galaxy hosts - 40 GCs (jKundu fc Whitmord l200l and 
we estimate the dynamical friction timescale before it 
disrupts and adds its stars and GCs to the main galaxy 
to be ~ 2 Gyr. 



5. CONCLUSIONS 

We have found radially-decreasing profiles of rotation 
and metallicity in the GC system of NGC 3115, and com- 
pared these to theoretical expectations. While the cen- 
tral regions were probably formed in violent, dissipative 
circumstances such as a gas-rich major merger, this event 
should have produced high or misaligned outer rotation. 
A more likely scenario for explaining the assembly of the 
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120101 ). Consequently, kinematical and chemical transi- 
tions such as those observed here should be ubiquitous. 
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